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Long-scale dynamic fluctuation phenomena in freely suspended films is analyzed. We consider 
isotropic films that, say, can be pulled from bulk smectic A liquid crystals. The key feature of such 
objects is possibility of bending deformations of the film. The bending (also known as flexular) 
mode turns out to be anomalously weakly attenuated. In the harmonic approximation there is no 
viscous-like damping of the bending mode, proportional to q 2 (q is the wave vector of the mode), 
since it is forbidden by the rotational symmetry. Therefore the bending mode is strongly affected 
by non-linear dynamic fluctuation effects. We calculate the dominant fluctuation contributions to 
the damping of the bending mode due to its coupling to the in-plane viscous mode, that restores the 
viscous-like q 2 damping of the bending mode. Our calculations are performed in the framework of 
the perturbation theory where the coupling of the modes is assumed to be small, then the bending 
mode damping is relatively weak. We discuss our results in the context of existing experiments and 
numeric simulations of the freely suspended films and propose possible experimental observations 
of our predictions. 


I. INTRODUCTION 

Thin flexible films and membranes (e.g., lipid bilay¬ 
ers) are curious but ubiquitous objects in the realm of 
condensed matter science. In the long-scale limit they 
can be considered as two-dimensional objects embedded 
into three-dimensional space. The main peculiarity of 
such objects is possibility of their bending distortions 
that have to be analyzed in addition to traditional long- 
scale deformations like density variations and in-plane 
motions. Here we turn our attention to freely suspended 
films that are not surrounded by a dense matter (unlike 
the biological membranes), that makes their bending os¬ 
cillations weakly attenuating and, as a result, sensitive 
to non-linear fluctuation effects. While basic thermal 
(and even quantum) fluctuation effects in various three- 
and two- dimensional condensed matter systems are well 
studied and known, the situation with freely suspended 
films is not completely recognized. In the paper we an¬ 
alyze liquid-like films that are isotropic thin sheets of 
matter. 

The very existence of freely suspended films, technol¬ 
ogy of their production and basic experimental and the¬ 
oretical studies are known from 70-ties of the previous 
century, and summarized in a number of review publica¬ 
tions and monographs (see, to quote some of them, M). 
Majority of experimental and theoretical investigations of 
physical properties of the free standing liquid-crystalline 
films and fluid membranes are devoted to thermodynamic 
and structural characteristics. In terms of correlation 
functions of long-scale degrees of freedom, it corresponds 
to measurements, numeric simulations or calculations of 
the static structure factor. Some works were performed 
also to study dynamic structure factor or time dependent 
correlation functions (see cited above literature and ref¬ 
erences therein), and only a very few results concerning 
non-linear dynamics and fluctuation effects are known 


M- However, in all those works the film was assumed 
either to have a finite (and not too small) surface ten¬ 
sion, or (for tensionless membrane) to be surrounded by 
viscous isotropic liquid. 

In a parallel world of crystalline or glassy freely sup¬ 
ported membranes investigations of the vibrational dy¬ 
namics also started long ago from the very influential 
works Bl. In these and some other papers (see also 
the review [3| ) in contrast to fluid-like films another situa¬ 
tion has been considered. Namely, solid tensionless mem¬ 
branes freely suspended in the vacuum (or in a dilute gas) 
were examined. These types of the investigations were 
resurrected recently (see e.g., fl3l - [T7| ) with a motivation 
(besides own fundamentally interesting physics) of vari¬ 
ous electro-opto-mechanical applications of the graphene 
films. However, fluctuation effects in solid and in liq¬ 
uid films are dramatically different. One can easily see 
the difference even in static properties. Indeed, as it is 
well known (compare results for solid films [a and fluid 
membranes |l0) or in [dj), thermal fluctuations make the 
freely suspended solid membrane more rigid, whereas for 
liquid membranes thermal fluctuations make the mem¬ 
brane more floppy. This qualitatively different behavior 
is accounted for the different in the both cases (solid and 
fluid membranes) sources of the anharmonic terms, re¬ 
sponsible for the bending fluctuations. In the former case 
the principal anharmonic terms come from the coupling 
of in-plane and out-of-plane (flexular) degrees of freedom, 
whereas in the latter case the bending softening occurs 
as a result of the intrinsic pure bending anharmonicity. 
Summarizing this brief overview of the state of art with 
film vibrational fluctuations, we conclude that there are 
some studies of fluid membranes surrounded by a viscous 
isotropic liquid and/or under tension, but the case of ten¬ 
sionless freely suspended in vacuum (or gas) membranes 
has been investigated only for solid membranes. The pur¬ 
pose of our work is to close this problem overlooked in 
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the previous works. 

An interest to fluctuation dynamic phenomena reflects 
both their practical importance and related fundamental 
science challenges. In this paper we extend the previous 
works, see Refs. , and also [9ij, to study dynamic fluc¬ 
tuation effects in the freely suspended in vacuum (or gas) 
tensionless liquid films. As experiments grow in resolu¬ 
tion and in sophistication, it is likely that further details 
of dynamics will be revealed whose interpretation will re¬ 
quire a reliable theoretical approach, we describe in our 
work. It is worth to note to the point that structurally 
similar to fluid films, polymerized or amorphous glassy 
membranes dynamically belong to the class of crystalline 
membranes (see, e.g., [16j). Indeed, in such membranes 
the in-plane overdamped viscous mode (whose coupling 
with the bending mode is the main source of fluctuation 
contributions into the mode damping) is replaced by the 
propagating transverse sound mode. 

Our paper is structured as follows. In the next section 
im we introduce main ingredients of our model, define 
physical quantities of interest, and derive basic relations 
for the quantities. In section mu we introduce the key 
tool for theoretical study of dynamic fluctuations - the 
effective action. Then in section HVl we give details for the 
perturbative treatment of the effective action functional, 
and formulate our main results in this work, namely the 
pure fluctuation bending mode damping and fluctuation 
contribution into the in-plane shear viscosity. In Con¬ 
clusion we present a summary of our results. Appendix 
contains some details of the calculations. 


II. BASIC RELATIONS 

We examine physical properties of freely suspended 
thin films. Very thin films (consisting of a few molec¬ 
ular layers) can be pulled from bulk smectic liquid crys¬ 
tal phases @, (H, [l9j] . It should be borne in mind that 
stability of such thin films in the freely suspended state 
is accounted for by the layered structure inherent in all 
smectics, i.e., by the presence of some internal interac¬ 
tion forcing molecules to stay in one layer. At increasing 
temperature, when this layered structure is destroyed, 
a freely suspended film becomes unstable with respect 
to arising holes induced by thermal fluctuations. There¬ 
fore the films exist in a restricted temperature interval. 
Another possibility to create very thin films is by using 
special dopes to conventional soap films the stability of 
which is supplied by surfactant monolayers on the water 
surface. Such films may be prepared with the thickness 
of several hundred Angstroms and investigated by opti¬ 
cal and X-ray methods (see e.g., references in and 

more recent papers |20l.l2ll]h 

There are different types of orientational ordering in 
the films. The films can be liquid (isotropic), such films 
are pulled from a bulk smectic-A phase. Films pulled 
from low-symmetry smectic phases possess an in-plane 
orientational order. The films pulled from bulk smectics- 


C have a “nematic” orientational order and ones pulled 
from bulk hexatic smectics possess a hexagonal orienta¬ 
tional order. Here we examine the simplest case of the 
isotropic films. 

We treat the films as two-dimensional systems, i.e., we 
assume that the thickness of the film is much smaller 
than the film lateral size. This confines our study to the 
scales, exceeding the film thickness. Let us stress that 
in studying freely suspended films there is no interaction 
with a substrate that is quite important for conventional 
two-dimensional systems formed on the surface of liquids 
or solids. Another characteristic peculiarity of such films 
is a possibility of their bending distortions. If the thick¬ 
ness of the film is sufficiently small fluctuation effects are 
noticeable on scales exceeding its thickness. The main 
goal of our work is to investigate the effects. 

We start our consideration with a derivation and anal¬ 
ysis of the film energy. If the surface tension a of the film 
is small one should take into account the contribution to 
the energy related to the film curvature in addition to 
the contribution from the film surface tension a. In the 
main approximation, the curvature contribution to the 
film energy can be written as 


H a = dS 


L2 


(R 1 1 + R 2 1 ) 2 + kR 1 1 R 2 


(1) 


by analogy with the bending energy of lipid bilayers (see 
the original paper pT)j ] and also its textbook version in 
!)■ Here R\ and i ?2 are local radii of the film and the 
coefficients k and R are called bending modules (or Hel- 
frich modules). The quantity R^R^ 1 is the Gaussian 
curvature of the film, and the combination R^ 1 + R 2 X is 
called its mean curvature. In accordance with the Gauss- 
Bonnet theorem, the last term in Eq. (JTJ) (with the coef¬ 
ficient R) is the topological invariant. Hence it does not 
play a role for fluctuations of the film shape that does not 
change its topology. That is why further (at examining 
fluctuations effects in the framework of the perturbation 
series) we ignore the topological contribution to the film 
energy. 

Comparing the curvature energy © with the energy 
related to the surfa ce tension, one finds the characteristic 
length l a = \Jk/<7. At scales larger than l a the surface 
tension dominates in the film energy. Dynamic properties 
of the film in the region of scales including fluctuation 
effects were examined in the work [7|. Note to the point 
that in an experimental setup, where the film area or the 
film thickness may relax to the equilibrium, the surface 
tension tends to zero since a = 0 is just the equilibrium 
condition. Such experimental setup is realized if the film 
is suspended without stretching or if the film is connected 
to a reservoir of molecules constituting the film. Then 
one expects that there exists the region of scales between 
the film thickness Ith and l a where the bending energy 
© dominates over the surface tension term. We examine 
further just the intermediate region of scales. The main 
goal of our investigation is to examine a role of dynamic 
fluctuation effects that play a relevant role in the long- 
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scale dynamic properties of the films. 

For scales larger than l a bending distortions of the film 
propagate, the mode can be termed as bending sound. Its 
velocity is Cb = (a / p ) 1 / 2 where p is the two-dimensional 
mass density of the film. The bending sound has the dis¬ 
persion law u) = Cbq, where ui is frequency and q is wave 
vector. For scales smaller than l a the bending distortions 
still propagate. However, the acoustic dispersion law in 
this region is substituted by the quadratic dispersion law 
u> = (k/ pY^ 2 q 2 ■ The g-dependence of the frequency in 
this case is similar to one for spin waves in ferromagnets. 

As in the case of the bending sound the bending mode 
with the dispersion law u> = (k/ p) r ^ 2 q 2 has an anoma¬ 
lously weak linear damping, because due to rotational 
symmetry viscous-like damping is forbidden. The rea¬ 
son is that if the bending viscosity is non-zero it would 
lead to an energy dissipation in a homogeneously rotat¬ 
ing film M, that is impossible. That is why the bare 
bending viscous coefficient is zero and the damping of the 
bending mode is determined mainly by non-linear fluc¬ 
tuation effects. In our case the main contribution into 
the damping comes from the interaction of the bending 
mode with the in-plane viscous mode (unlike crystalline 
membranes [ljf], or liquid crystalline films embedded in 
a viscous liquid 

In the region of negligible surface tension, for scales 
between Z t h and l a , thermal bending fluctuations of the 
film shape lead to a logarithmic renormalization of the 
modules k and R. First an attempt to calculate the renor¬ 
malization of the module n was taken by Helfrich [H |. 
and later by Forster (24j . The correct renormalization 
group (RG) equation for the bending module k in the 
one-loop approximation was derived in papers by Peliti 
and Leibler [2l|, Kleinert [2(| and Polyakov {27(. The 
RG-equation for R in the same approximation was found 
by Kleinert [28j. The RG-equation for the bending mod¬ 
ule k is 


dn/dt, = — 3T/(4tt), (2) 


monograph i) concerning the fluctuation damping of 
the hydrodynamic modes in freely suspended films for 
the region of scales larger than l a . The bending sound 
damping comes mainly from thermal fluctuations, the 
damping can be estimated as Tq 3 /(pcb). There is also a 
fluctuation contribution to the damping of the in-plane 
sound that can be estimated as 

Z r ! ) 1/3 ' <3) 

However, because of a small numerical factor (on the or¬ 
der of 10 -2 ) in front of the quantity (0 in the sound 
damping one should take into account the conventional 
damping ~ pq 2 / p in addition to the fluctuation contribu¬ 
tion. Analogously, one should take into account the con¬ 
ventional damping pq 2 /p of the transverse (to the wave 
vector q) velocity perturbations in addition to the fluc¬ 
tuation contribution f3]) (with a small numerical factor 
on the order of 10~ 2 ). 

Note that the specific entropy s/p (where s is the two- 
dimensional entropy density, s dx dy is the entropy of the 
film element) is hardly excited by the fluctuations we are 
investigating. That is why below the quantity is assumed 
to be equal to its equilibrium value. Therefore further on 
all thermodynamic derivatives are assumed to be taken 
at constant specific entropy s/p. 

We assume that in equilibrium the film is parallel to 
the X — Y plane and its bending distortions are charac¬ 
terized by the displacement h of the film along the in¬ 
direction . The displacement h is treated as a function of 
x , y. Therefore the film shape is determined in the Monge 
representation as z = h[x,y). Then the area element of 
the film is expressed as dS = g 1 ! 2 dx dy, where 

g — 1 + ( d x h ) 2 + (d y h) 2 . (4) 

The quantity g plays a role of the determinant of the 
metric tensor of the film. The unit vector perpendicular 
to the film has the following components 


where £ is logarithm of the scale, where the modulus is 
determined, and T is temperature (measured in energy 
units). The approach implies that T/k <C 1. 

We see that the dimensionless parameter T/n charac¬ 
terizes intensity of the thermal bending fluctuations. For 
lipid bilayers, the quantity is usually on the order of 10 -2 . 
The smallness is related to the fact that the bilayer thick¬ 
ness exceeds the atomic length. Therefore one expects 
that for the freely suspended films the ratio T/k is even 
smaller since the bending module k is roughly propor¬ 
tional to the third power of the film thickness (as it fol¬ 
lows from the classical theory of elastic shells). Further, 
we treat the ratio T/k as a small parameter. There is 
another dimensionless parameter characterizing the film, 
np/p 2 , where p is the in-plane shear viscosity coefficient 
of the film. Like it is in the case for conventional bulk 
nematics, we expect that the parameter np/p 2 is small 
as well. 

Let us recall main results of the work (7J (see also the 


k=9 1/2 {-d x h,-d y h, 1). (5) 


The starting point of our analysis is the energy of the 
film. Further, we treat all variables characterizing the 
film as functions of x , y. Then the first term in the energy 
© is rewritten as 


'H s = dxdyg 1/2 ^ d a (d a h/g 1/2 ) 


( 6 ) 


Here and below Greek subscripts designate components 
of the vectors along the X , Waxes. In addition to the 
bending energy ©, we consider also the film energy re¬ 
lated to the film compressibility 

H n = Jdxdyg 1 ' 2 ^-^-?) , (7) 

where B is the compressibility module. Here p 0 is the “in¬ 
trinsic” equilibrium mass density of the film and pog 1 ^ 2 is 
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the equilibrium mass density in projection to the X — Y 
plane. Note that the expressions (Em are formally valid 
for arbitrary variations of the film shape, the only re¬ 
striction is absence of over-folds, that is uniqueness of 
the function h(x,y). The contributions (16171) have to be 
supplemented by the kinetic energy of the film that is 
written as 

r I 2 -)- j 2 -|- j 2 

H kin = J dxdy Jx 3 * Jz , (8) 


Here g and £ are two-dimensional analogs of the three- 
dimensional first (shear) and second (bulk) viscosity co¬ 
efficients. 

In the linear approximation (for the region of scales 
under consideration, smaller than l a ) one finds the fol¬ 
lowing hydrodynamic (i.e., long-scale) modes of the film: 

• the in-plane propagating sound mode with the lin¬ 
ear dispersion law co = ciq ; 

• the overdamped in-plane thermodiffusion mode; 


where j is the film momentum density. Its definition 
implies that j dx dy is the momentum of the film clement. 
Then the film velocity is v = j/p, as usual. 

The film surface tension a is a variation of the film 
energy over the film area. Taking into account that the 
area element is g 1 ^ 2 dx dy and varying the expression GD 
over g 1 / 2 , one obtains 


a = — B 




( 9 ) 


that is an expression of the first order in the film den¬ 
sity deviations from its equilibrium value. If the film is 
stretched then its mass density p deviates from the equi¬ 
librium value pag 1 / 2 and the surface tension is non-zero. 
A non-zero surface tension is produced by fluctuations of 
p. Note that just the module B determines the in-plane 
sound velocity eg. c 2 = —da / dp = B / p. 

The dynamic non-linear equations for the film displace¬ 
ment h and for the mass density p are 

pdth — jz jadah, (10) 

dtp = - daja • (11) 


The equation m is the kinematic condition implying 
that the film moves with the velocity v = j /p, and the 
equation m is the mass conservation law. Therefore the 
equations (110I11D are exact that is there are no dissipative 
corrections to the equations. 

The equations (1101111) should be supplemented by the 
equation for the momentum density j£ of the film that is 
written as the momentum conservation law 


dtji — d a ^'Vaji T g ^ T,a 9 ^ giafimdpVm^ , (12) 

where Latin subscripts run over x 7 y 7 z. All quantities in 
the equations are, as above, assumed to be functions of 
x,y. The expression for the stress tensor Tik was found 
in the work [8] (see also the book i) , it is written as 


%k 


a+^-{Xlf 


5i + K{X7l)d£h 


nhd^Xl. (13) 


Here 5^. = <5^ — IjJk stands for the projector to the film, 
the unit vector l has the components ([5]), Vi = d a l a7 and 
d£ = dkada- The viscosity tensor can be written as 


Viklm = {V~ 0 5a S L + V{Sti S km + S L S ki)- (I 4 ) 


• the overdamped in-plane viscous mode; 

• the bending (flexular) propagating mode with the 
dispersion law u> = yj k/ p q 2 . 

The specific entropy dynamics (thermodiffusion mode) 
can be separated from the other modes. Then, in what 
follows we assume, that the specific entropy mode is not 
excited, and therefore one can safely forget about the 
thermodiffusion. Note that in the linear approximation 
the in-plane sound mode has the standard viscous damp¬ 
ing ~ ( g/p)q 2 , whereas for the bending mode such kind 
of viscous damping is forbidden bv the rotational sym¬ 
metry of the film (for details see [9J). Only higher order 
over the gradients (e.g., proportional to g 4 ) dissipative 
terms are not forbidden, and such terms will produce a 
non-zero damping of the bending mode. However, as we 
demonstrate below, thermal fluctuations produce a larger 
contribution to the bending mode damping, that is pro¬ 
portional to q 2 . Calculation of the damping is one of the 
main goals of our work. 

In the hydrodynamic regime (that is for scales, larger 
than the film thickness) the in-plane sound is harder, i.e., 
its frequency co for a given wave vector q is much higher 
than the characteristic frequencies (inverse time scales) 
for the bending and for the in-plane viscous modes. 
Therefore the degrees of freedom related to the sound 
(density fluctuations and the longitudinal to the wave 
vector component of j a ) can be effectively excluded from 
the consideration. Then for the first step of the analysis 
we stay with the closed description of the bending and 
of the shear viscous modes. We analyze their non-linear 
interaction leading, particularly, to fluctuation contribu¬ 
tions to damping of the modes. 

Particularly, we establish, that the non-linear fluctua¬ 
tion effects produce the viscous-like damping of the bend¬ 
ing mode that leads to the following dispersion law 

co = \Jn/p q 2 - iuq 2 /2, (15) 


where v is the fluctuation bending viscosity. The second 
term in the right-hand side of Eq. (1151) gives the bending 
mode damping related to scattering of the bending mode 
on fluctuations of the viscous in-plane mode. That is why 
v is proportional to the temperature T: 


In 

7TK,p 


(16) 


A derivation of the expression (1161) is given in the next 
sections. 
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III. EFFECTIVE ACTION 

We proceed to calculating the bending mode damp¬ 
ing 1151161) and the corrections to the shear viscosity q 
caused by thermal fluctuations. The calculations can be 
done in the framework of the diagrammatic technique 
first developed by Wyld [29j] in the framework of hydro- 
dynamic turbulence. Then the technique was generalized 
for a wide class of systems by Martin, Siggia and Rose 
[30| . Note that the technique is the classical limit of the 
Keldysh diagrammatic technique 0] (see also the text¬ 
book B). We use the version proposed by Janssen jH], 
where correlation functions of the fluctuating fields are 
written as path (functional) integrals over the fields with 
the weight exp (*/) where / is an effective action con¬ 
structed in accordance with the dynamical equations of 
the system (see also the book 0). 

The effective action / for the problem we are investi¬ 
gating is the sum / = / reac + /diss of the reactive and of 
the dissipative contributions, 

-^reac — J' dt dx dy [p(pdth j z T j a d a h) 

T Pi^tji T Pi^cx.{Pcx.ji) T Pi^a (v^Tia:)] ’ (1^) 

-^diss — j ^ ^ daPiTliaPm [dfjVm + ’iTdpPm] , (18) 

where p and p are auxiliary fields. Here we omitted the 
specific entropy fluctuations. The action I is constructed 
in accordance with the dynamical equations m and 
ed , and the relation (fill) is implied. The last term in 
Eq. 1181) comes from viscous Langevin (random) forces. 
That is why the term is proportional to the temperature 
T. Note that the action 1171181) is formally exact in h. 

Note that the pair correlation functions, like (hp) (here 
and henceforth angular brackets mean averaging that 
is equivalent to a functional integral with the weight 
exp(*/)) determine response of the system to an exter¬ 
nal excitation (force). The structure of the action (fl8l) 
guarantees validity of the fluctuation-dissipation theo¬ 
rem. The theorem leads to the relation between the cor¬ 
relation functions like ( hh ) and {hp). Note also, that the 
pair correlation functions like (pp) are zero. 

Both, the bending and the in-plane viscous modes have 
the bare dispersion laws w oc q 2 . Therefore in the long- 
scale limit (small wave vectors) in the main approxima¬ 
tion one can neglect the time derivative in the mass con¬ 
servation law HID to obtain the condition d a j a =0. It 
is analogous to the incompressibility condition for the 
Navier-Stokes equation. Next, an integration over the 
density field p (that is the variable related to the in-plane 
sound) at calculating the correlation functions can be 
performed in the saddle-point approximation. By other 
words, the action I can be substituted by its saddle-point 
value that can be found by equating to zero the variation 
of the action 1171) over p. The main p-dependence of the 
action is related to the contribution to the stress tensor 
Tia, proportional to the surface tension a, see the expres¬ 


sions 191131). Calculating: the variation of the action 03 
over p and equating the variation to zero, one obtains the 
following condition 

Sai d aPi = 0. (19) 

Particularly, the condition (flT)l) leads to the conclusion 
that the term with the surface tension and the term with 
the difference rj — £) in the viscosity tensor ED do not 
enter the effective action for the interacting bending and 
viscous modes. 

In this approximation, the surface tension a is not a dy¬ 
namic variable, like pressure in the Navier-Stokes equa¬ 
tion. The quantity a is passively followed to the bending 
and the viscous fluctuations. Thus, variations of the sur¬ 
face tension are relatively weak. Therefore, to determine 
the mass density, the surface tension a can be put to 
zero. Then, we obtain from Eq. © that 

p = p 0V / l + (Vh) 2 , (20) 


where po is the “internal” mass density of the film. The 
expression © leads to the following renormalization law 


d In po T 

dt ; AlTK ’ 


( 21 ) 


where £ is logarithm of the current scale. 

One should choose variables describing the bending 
and the viscous modes. We take the physical variables 
h,j z ,j a . The last one satisfies the condition d a j a = 0, 
therefore we deal with three scalar fields. As to the aux¬ 
iliary fields, we choose p, p z and the transverse (to the 
wave vector) component of p a , satisfying the condi¬ 
tion d a Pa = 0. 

In the framework of the perturbation theory one starts 
from the “bare” correlation functions that can be easily 
obtained by Gaussian integration if we keep solely the 
second-order term in the effective action I. Explicit ex¬ 
pressions for the correlation functions of the viscous mode 
are given by the expressions 


(j a (t,r)jp( 0,0)) = J 

^ __ 2Tr)(q 2 6 a p - q a q 0 ) 

" oc/3 o i 2 4/2 1 

or + rj z q^/p z 

Oa(t,r)p£(0,0)> = j 

n _ q a qp/q 2 

^*0.(3 2 / 

uj + irjq z /p 


( 22 ) 


(23) 


However, to analyze fluctuation effects one should use 
“dressed” correlation functions of h 7 j z , i.e., with in¬ 
cluded fluctuation contributions. Then the pair corre¬ 
lation function of the displacement h is written as 


(h(t,r)h( 0,0)) = 
Fhh = 


dui d 2 q . . 

exp(—zwf + iqr)F hh , 

(270 

2n (w,q) 


[w 2 - nq 4 /p+ £(w)][w 2 - nq A /p + £(—w)] 


,( 24 ) 
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where II and E are “polarization” and “self-energy” func¬ 
tions, in the terminology borrowed from the quantum 
field theory. The quantities II and E have to be calcu¬ 
lated in the framework of the perturbation theory. Due 
to the fluctuation-dissipation theorem 

Im E(w, q) = -Im E(-w, q) = y^II(w, q). (25) 

One finds that the main contributions to II and E are 
proportional to q 2 . The real part of E reproduces the 
renormalization of n and p 0 that is irrelevant for us due 
to the inequality T/k <C 1. That is why further we take 
into account only the imaginary part of E: ImE = vq 2 , 
II = (T / p)vq 2 . In the approximation we obtain from Eq. 

(El 


(hh) = J 


dco d 2 q , . 

- — exp (~iut + iqr)F hh , 


(2tt) 

Fhh = 


2 Tp 1 vq 2 


(ui 2 — nq 4 / p) 2 + v 2 q 4 u> 2 
The expressions correspond to the dispersion law (fl5l) . 


(26) 


IV. PERTURBATION THEORY 


In this section we calculate the non-linear fluctuation 
effects in the framework of the perturbation theory. For 
the purpose we expand the effective action I. We use 
the one-loop approximation that is one should expand 
the effective action up to the fourth order. Results of the 
calculations in the framework of the perturbation theory 
can be presented by Feynman diagrams, where lines cor¬ 
respond to the pair correlation functions (I22I23I26I) and 
vertices are determined by the third and fourth order 
terms in the effective action. Note that the expansion 
of the contributions to the action (11711811 has to be per¬ 
formed in h since the terms are quadratic in the other 
fields. One should expand in h the relation (fl9l) as well. 
The expansion determines the longitudinal (to the wave 
vector) component of p a . We find 

d a p a = —d a hd a p z + d a hdphd a pp, (27) 

with the second-order accuracy in h. 

Now we proceed to calculate the fluctuation contri¬ 
bution to the polarization function II entering the ex¬ 
pression EH- Because we assume rj 2 up the leading 
interaction vertices are determined by solely the dissipa¬ 
tive effective action m- To find II one should select the 
terms with the field p z . All relevant third order terms 
containing p z are 


/third = 


dt dx dy | ^ dphd a p z (dpj a + d a jp) 
+2iTrjdphd a p z (dpPa + daPp')^. (28) 


The fourth order term needed for us is 

/fourth = J dtdx dy 1 2iT rjd jt ,hd v hd tl p z d v p z 

-AiTi]d fJ ,p z d^h^^-(dphd ol p z ) 

+iTr)[d/ 3 hd a h + 5cq3(V/i) 2 ]<9 a p z <9/3p 2 j, (29) 

it is quadratic in p z . Note that the action (051) is non-local 
(it contains the term proportional to 1/V 2 ), as one could 
expect because this action is derived by the integration- 
out of the full effective action over the acoustic degrees 
of freedom. 



There are some one-loop contributions to the quantity 
II, that can be represented by the Feynman diagrams, see 
Eqs. (130113211 . The cross in Eq. EHl) represents the fourth- 
order vertex determined by the action 051) and the wavy 
line stands for the pair correlation function (12611 . The 
third-order vertices in Eqs. EM are determined by 
the action (051) . The solid line there represents the pair 
correlation function (051) and the combined solid-dashed 
line represents the pair correlation function 031) . The 
diagrams enable one to write an explicit expression for 

n. 

The structure of the vertices leads to the conclusion 
that n(fc) oc k 2 , where k is wave vector. Substituting 
II = (T / p)vk 2 one finds after lengthy but straightforward 
mathematics 

vnp _ k a kp f d 2 q (q a + k a )(qp + kp) 

Tp k 2 J (27r) 2 \q + fc| 4 

_ 4 MV f d 2 q (q^ + ky)(qp + kp) q a qp 

k 2 J (27r) 2 |q + k\ 4 q 2 

kpk u f d 2 q (q a + fc a )(g M + fc M ) 1 
k 2 J (27t) 2 \q + k\ 4 q 4 

{daqv{q QpQfi) T Q/3Qv(.Q &afj. daq^,) 

Tqaq^i[q dftv qpqu ) T qpq^i(q dav qa.qv )}■ (33) 

Simple inspection of this expression shows that all the 
integrals converge in the both limits, at large q and at 
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small q + k. Therefore the main contribution into the 
integral comes from the region q ~ k and there are no 
ultraviolet contributions (particularly, logarithms) to the 
quantity v, in accordance with the symmetry arguments. 
The dimension analysis leads to the conclusion that the 
quantity ffTTTTI) is independent of k. Calculations (collected 
in Appendix) give the final answer (1161) . 

Out of these calculations will come a bit unexpected re¬ 
sult, that forbidden by rotational symmetry viscous-like 
damping of the bending mode, becomes perfectly legiti¬ 
mate driven by thermal fluctuations. Note to the point, 
that the situation where the fluctuation damping (de¬ 
termined by fluctuations of the same scale as the wave¬ 
length) has some scaling forbidden for the bare damping 
by symmetry is characteristic also of the spin waves in 
two-dimensional ferromagnets [37]. 

Let us estimate a fluctuation correction to the shear 
viscosity coefficient q. As previously, we consider only 
the interaction terms produced by the dissipative part of 
the effective action m, that is justified by the inequality 
up <C q 2 ■ One obtains the principal interaction term 


Ant = J dtdxdy [{r]/p)dphd a j z (dpp% + d a pf) 

+2iTr]dphd a p z (d l 3p t a + d a pf)}. (34) 

The one-loop contribution to q is determined by the same 
diagrams (1311321) where interaction triple vertices are de¬ 
termined now by the action (15311 . Calculations lead to 
the following estimation of the fluctuation correction 


Arj ~ 



(35) 


Now we can accurately formulate the applicability con¬ 
ditions of our approach. Our procedure is correct pro¬ 
vided Ary -C q that is if 


e 


/ q 
K yfKp 


< 1 . 


(36) 


Note that the condition (1551) implies that the damping of 
the bending mode is much less than its frequency since 
the inequality (1551) is equivalent to the condition v -C 
\/n/p , see Eq. (TTSl) . 

As a note of caution we should say also, that above 
we neglected the self-interaction of the viscous mode that 
leads to the logarithmic corrections to the shear viscosity 
coefficient ry that can be estimated as [HI, [3j| 


Ary ~ pT/q, 


up to a logarithmic factor. Comparing the estimation 
with Eq. (1551) we conclude that the effect we calculated 
above is stronger due to the inequality ry 2 up. 

It is interesting that just the parameter (1551) is the 
small parameter justifying the perturbation expansion in 
h. If the parameter is not small then we are beyond 
the applicability of the perturbation theory, in the re¬ 
gion of strong interaction. Then nothing can be done 


in the framework of the perturbation theory. Based on 
heuristic hand-waving arguments related to the fluctu¬ 
ation enhancement of the viscosity coefficient we could 
only speculate that probably, the behavior in that region 
corresponds to a glass state of the film since an essential 
viscosity growth is expected in the case. 


V. CONCLUSION 

In conclusion, in this paper we describe dynamic fluc¬ 
tuation phenomena in freely suspended films. We re¬ 
strain ourselves to liquid-like (isotropic) films in the 
overlooked in the previous works conditions (tensionless 
membrane freely suspended in the gas or vacuum). The 
bending (flexular) mode of such film turns out to be 
soft and weakly attenuated, therefore yields to strong 
dynamic fluctuation effects. In the harmonic approxi¬ 
mation, for the bending mode there is no viscous-like, 
proportional to q 2 (q is the wave vector of the mode) 
attenuation, and only much smaller super-viscous atten¬ 
uation oc q 4 is not forbidden by the rotational symmetry. 
We calculate the dominant fluctuation contributions to 
the damping of the bending mode due to its coupling 
to the viscous (non-propagating) in-plane mode. The 
fluctuation damping restores the viscous-like q 2 attenu¬ 
ation law of the bending mode. The damping is weak 
due to smallness of the dimensionless coupling constant 
obligatory within our perturbation approach. What is 
fascinating about our results is that they not only con¬ 
tribute to understanding of many dynamic properties im¬ 
portant for nano-technological opto-mechanical applica¬ 
tions of the films, the results can be confronted with re¬ 
cent experimental data [3(1 on fluctuation enhancement 
of the membrane viscosity. Another way to check exper¬ 
imentally our theory predictions is to perform inelastic 
(dynamic) light scattering experiments (like in the works 
I 2 QU 2 H) to measure the scattered signal (intensity) line- 
width. 

What makes our calculations involved and non-trivial, 
is the fact that there are two essential dimensionless pa¬ 
rameters in the theory, namely, T/ k <Cl and q 2 / (up) 

1, and physics depends crucially on how they interplay. 
The small parameter of our perturbation theory is e, see 
Eq. (1551) . One can try to extract the value of e from liter¬ 
ature data on material properties of the lipid membranes 
and the liquid crystals that are rather dispersed, see, e.g., 
I3H . There are some materials where (at room temper¬ 
atures!) e < 1, and also those with e > 1, but the both 
inequalities are not too strong. One also should keep in 
mind that there is regular temperature dependence (a 
sort of Arrhenius law) of the bare material parameters 
[3§| . If we stretch all essential material parameters to 
their utmost (but still not unrealistic) values, the ratio 
could be either e < 1, or e > 1. If the ratio is not small 
then we go beyond the applicability of the perturbation 
theory, into the region of strong interaction. This prob¬ 
lem deserves a separate investigations. In this regime 





of very strong bending fluctuations one can expect the 
membrane to be broken via non-perturbative processes, 
like pore formation. Another theoretically tempting pos¬ 
sibility would be vitrification of the membrane. 

We believe that our results are also a step forward 
in establishing connections between different regimes of 
the collective membrane dynamics. In own turn think¬ 
ing about membrane dynamics from a pure physics-based 
ground can bring new insight on some relevant character¬ 
istics of biophysical cells. We expect fluctuation dynam¬ 


ics of freely suspended membranes to be a fascinating 
and fruitful field of investigations for biologists, applied 
scientists, and experimental and theoretical physicists. 
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Appendix A 

To calculate the integral (1531) we substitute there k = 
(1, 0) to obtain 

"up _ f dip dqq _1_ 

Tq J (27r) 2 (1 + q 2 — 2q cosip) 2 

5 1 

q 2 cos(2<^) — 2 q cos ip + — — 2 cos(2^>) + - cos(4^>) . 
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The integral over the angle ip can be found using the 
following relations: 

1 f 2n dp cos (mip) 

2ir J 0 (1 + q 2 - 2q cos ip) 2 

(to + 1 )q m + (1 — m)q 2+m 
= ’ 

if 0 < q < 1, to > 0, and 

1 r 2w dip cos (mip) 

2i r J 0 (1 + q 2 - 2q cos ip) 2 

(to + 1 )q 2 ~ m + (1 — m)q~ m 
= ’ 


if q > 1, to > 0. Thus the integral over q is divided 
into two intervals: from 0 to 1 and from 1 to oo. The 
integration can be easily performed and one finds 


+ 


UKp 

W = J 0 

f°° dq q 


1 dq q 


1 


2-7T (1 — q 2 ) 3 2 
1 3 


(1 — h) 


2\3 


2tt (q 2 — l) 3 2q A 


(q 2 -iy = -. 


(A2) 


Thus we arrive at the result (USD- 








